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Photoionization by attosecond (as) extreme ultraviolet (xuv) pulses into the laser-dressed con- 
tinuum of the ionized atom is commonly described in strong- field approximation (SFA), neglecting 
the Coulomb interaction between the emitted photoelectron (PE) and residual ion. By solving the 
time-dependent Schodinger equation (TDSE), we identify a temporal shift St in streaked PE spec- 
tra, which becomes significant at small PE energies. Within an eikonal approximation, we trace 
this shift to the combined action of Coulomb and laser forces on the released PE, suggesting the 
experimental and theoretical scrutiny of their coupling in streaked PE spectra. The initial state 
polarization effect by the laser pulse on the xuv streaked spectrum is also examined. 

PACS numbers: 42.65.Re, 79.60.-i, 42.50.Hz, 



I. INTRODUCTION 

Using an xuv attosecond pulse to photocmit electrons 
from gaseous or solid targets into the electric field of a 
synchronized delayed femtosecond (fs) infrared (ir) laser 
pulse provides a powerful tool for investigating ultrafast 
electron dynamics by recording ir-laser-streaked xuv PE 
spectra [l(. For strong laser fields and sufficiently fast 
PEs, streaked PE spectra are conveniently described in 
SFA Q, i.e., by ignoring the interaction of the residual 
ion with the released PE. In this case, subject only to 
the ir laser electric field, the propagation of PEs can 
be described in terms of analytically known "Volkov" 
states Q. This leads to a delay-dependent energy shift 
S e coe( t ) = -kA L (-r) of the center-of-energy (COE) 
Ecoe = ojx — \eb\ in the PE spectrum, where u>x is the 
xuv-pulse central frequency, Eb the binding energy in the 
initial bound state, k the PE asymptotic momentum, and 
At, the vector potential of the ir pulse. We use atomic 
units except where stated otherwise and define the delay 
r between the centers of the xuv and ir pulse as positive 
if the xuv pulse precedes the ir pulse. 

The interpretation of sub-fs temporal shifts in streaked 
PE spectra is a matter of current debate. For exam- 
ple, the recent measurement Q of a relative delay of 
« 1 10 ± 70 as between the ir-streaked xuv photoemission 
from localized 4f core levels and dclocalized conduction- 
band (CB) states of a W(110) surface was understood in 
the original ref. [H and a subsequent theoretical work [j| 
as the difference 8tcB-if = tcB — tif between the 
arrival times of CB and 4f core PEs at the surface. 
This interpretation is based on the assumption that the 
ir pulse does not penetrate the surface, such that CB 
and 4f electrons that are released at the same time — r 
by the absorption of an xuv photon get streaked only 
upon arrival at the surface, producing the COE shifts 
Ecoe = —kA L (—T + tcB) and -kA L {—r + Uf), respec- 
tively, in the PE spectra. According to this two-step 
explanation (photorelease followed by streaking) , the to- 
tal emission probability Ptot = J dEP(E, t) from a given 
initial state would not depend on r. 



In contrast, our analysis of experimental streaked pho- 
toemission data for a tungsten [j] and rhenium surface @ 
indicates that Pcb(t) oscillates with Al(t), with an am- 
plitude of £3 10% of the average value. Furthermore, the 
continuity of the wavefunction and its derivative at the 
surface, implies that an intense fs ir pulse affects the 
PEs inside the solid, even if the ir electric field were 
prohibited from penetrating the surface Q. We have 
shown that this observed temporal shift can be repro- 
duced within the SFA and interpreted it as an interfer- 
ence effect in the emission from different lattice sites 0], 
observing that the SFA cannot account for relative tem- 
poral shifts in the emission from different levels of iso- 
lated atoms. A classical transport simulation including 
the effect of (in)elastic collisions of released PEs with 
tungsten cores on the propagation of PE inside the solid 
leads to 5t C B—if = 33 as Q. Thus, different models [3- 
H, H strongly deviate with regard to the assumed at- 
tenuation of the streaking ir electric field Et, (t) inside the 
solid, ranging from no penetration into the surface Q to 
penetration depths of 30 0] and 85 layer spacings 
or larger than the electron mean free path Q. The de- 
tailed modeling of the (relative) delay in the photoemis- 
sion from metal surfaces is further complicated by the 
complex band structure and the ensuing difficulty in as- 
signing a group velocity to the motion of PE wave packets 
inside the dispersive conduction band surface charge 
accumulation, and the general concern that static band- 
structure calculations and the assumption of an instan- 
taneous plasmon response (i.e., static image charge in- 
teractions) are invalid at the as time scale. These shifts 
may be of particular importance in the interpretation of 
streaking spectra for complex targets, such as metals, and 
emphasize the need for more detailed studies of streaked 
photoemission spectra. 

It is of fundamental importance to first understand all 
contributions to this temporal shift for simple systems. 
In this work, we focus on the effect of simultaneous ir 
laser pulse and Coulomb interactions on streaked pho- 
toemission spectra from the prototypical ground state of 
a one-dimensional hydrogen atom. This Coulomb-laser 
coupling effect was first investigated by Kroll and Wat- 
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FIG. 1: (Color online) Streaked photoemission from ID model 
hydrogen atoms. TDSE calculations for xuv pulses with (a) 
fkdx = 90 and (b) 25 eV. (c) Corresponding centers-of-energy 
SEcoe{t) for ftwx = 90 (solid line) and 25 eV (dashed line). 
To facilitate the identification of the relative temporal shifts 
St, SEcoe{t, hwx = 90 eV) is normalized to the ftojx = 25 eV 
result, (d) 5t and (e) oscillation amplitude relative to the SFA 
for TDSE (full line) and eikonal approximation (dashed line) 
calculations. 

son [l(| in their study of laser-assisted atomic scattering. 
It also affects the spectra of high harmonic generation, 
multiphoton ionization, and laser-assisted xuv photoion- 
ization. For example, the Coulomb interaction causes the 
xuv streaked PE spectra to be right-left asymmetric [ll[ . 
In RABITT measurements (Reconstruction of Attosec- 
ond Beating by Interference of Two photon Transition), 
these simultaneous ir laser pulse and Coulomb interac- 
tions induce a so-called atomic phase which shifts side- 
band intensities as a function of the delay between the 
xuv pulse train and the ir pulse fl2l - [T5j . 

In this work, we demonstrate, numerically and ana- 
lytically, how the coupling of the ir laser pulse and the 
final-state Coulomb interaction of the PE with the resid- 
ual ion gives rise to a significant temporal shift St in the 
COE of streaked PE spectra with respect to those ap- 
proximated in SFA. As we will show, inclusion of this 
Coulomb-laser (CL) coupling alters both, the amplitude 
and phase of the COEs in streaked PE spectra, leading 
to a COE shift SE^q E (t) = -KA l (t - St) with an os- 
cillation amplitude K > k. Thus, St and the streaking 
amplitude ratio K/k i) help to reveal details of the PE 
dynamics including the combined interaction of Coulomb 
and laser forces and ii) converge to their SFA limits, 
and 1, respectively, at sufficiently large PE energies. 

We numerically solve the TDSE, including the 
electron-proton interaction, and then compare our re- 
sults for P(E, t) with SFA and eikonal approximation 
(EA) calculations for a large range of PE kinetic ener- 
gies. Strongly dependent on wx, we find temporal shifts 
St of more than 50 as. Fig.s Q] (a) and (b) show the 
TDSE ir-streaked PE spectra P{E, t) for gaussian xuv 
pulses of length tx = 300 as with fwjx = 90 and 25 eV, 
respectively. The ir pulse is also assumed to be a gaus- 
sian and has a peak intensity of II = 2 x 10 12 W/cm 2 , 



a carrier frequency u>l = 1-6 cV/h, and a pulse length 
t l = 5 fs. These spectra are shifted by St = 60 as, 
which becomes apparent in the corresponding COE shifts 
SEcoe(t) in Fig. Q] (c). The solid curves in Fig.s [1] (d) 
and (e) show St and the ratio Kj k of streaking-oscillation 
amplitudes for a large range of ujx- Within an EA ap- 
proach [ll|j [l6|, [I?]], we can trace (details will be given 
further below) this wx-dependent temporal shift and the 
oscillation amplitude enhancement to the CL coupling in 
the PE final state (dashed curves in Fig.s[T](d) and (e)). 
St and K/k, for the TDSE and EA calculations, con- 
verge at large ujx to their respective SFA limits St = 
and K/k = 1 due to the diminishing influence of the 
residual ion's Coulomb force at increasing PE energies. 
Thus, St and K/k — 1, are measures for the combined 
action of the Coulomb and laser force on the PE relative 
to the action of the ir laser force alone. Since St < 0, the 
attractive Coulomb force does not delay the PE emission, 
as one might intuitively expect. We also note that K > k 
reveals a Coulomb-enhancement effects that is reminis- 
cent of the Coulomb potential's infinite range leading to 
well-understood "Coulomb-Cusps" in energy-differential 
collision-induced PE spectra |18j . 

This article is organized as follows. In Sec. HJ we 
present numerical results based on the TDSE. In Sec. 
IIII1 we adopt an EA to take into account the simultane- 
ous ir and Coulomb interactions of the PEs, and com- 
pare our EA and TDSE results. In Sec. IIV( wc examine 
the effect of the polarization of the initial state on the 
streaked spectrum. We conclude in Sec. V. In the ap- 
pendix, we show that, within the eikonal approximation, 
the obtained atomic phase in RABITT is identical to the 
relative temporal shift in the streaked PE spectrum. 



II. TIME-DEPENDENT SCHRODINGER 
EQUATION FOR STREAKING 

The exact wavefunction of the one-dimensional model 
atom interacting with the ir and xuv pulse is determined 
by the TDSE (in the length gauge) 



2 dx 



, U{x) + V{x,t) 



*(x,t), (I) 



where U(x) is the Coulomb potential, V(x,t) = 
x[EL{t) + Exit)] the interaction with the ir and xuv 
pulse, and E L ( X ) the electric field of the ir (xuv) pulse. 
Assuming single-photon ionization in a sufficiently weak 
xuv pulse, and after splitting the exact wave function 
for the atom in the combined xuv and ir electric fields 
according to *f>{x,t) = ip g (x,t) + SxjAx,t), ([T]) can be re- 
placed by two coupled equations [l9( 



l-glpg{x,t) 



i-Q t S^{x,t) 



1Pg(x,t), (2) 

5ip{x, t) 



3 



+ XE X (t + T)lPg(x,t). 



(3) 



Equation © determines the evolution (polarization) of 
the initial state in the ir field and ([3]) the generation of 
PE wave packets by the xuv pulse and their evolution in 
the ir field. The electric fields £x(X) = —dA L i X ){t)/dt of 
the ir (xuv) pulses arc derived from the vector potentials 

A HX )(t) = A L{xho cos(LO L{x} t)e- 2los2t2 ^K Since 
EL(x)(t — > ±oo) = 0, equations (0 and are subject 
to the initial conditions ip g (x,t — > — oo) = tp(x)e~ leBt 
and Stp(x,t — > — oo) = 0. The ground-state initial wave 
function tp(x) and energy eb are obtained from 



e B il>(x) = 



Id 2 



ij){x). 



(4) 



We solve ©-(H) numerically by wave-packet propagation 
for times \t\ < 2.5tx with a step size At = 0.2 on a spatial 
grid with \x\ < 2000 and spacing Ax = 0.25. Assuming 
free-electron dispersion, E = ifc 2 , we calculate the ir- 
assistcd xuv photocmission probability 



P(E, 



8^p{k, r, t ■ 



(5) 



and the corresponding COE 
1 



E, 



COE 



(r) 



dk 



k8tp(k,r,t ^oo) /P tot (r), (6) 



where dip(k, r, i) is the Fourier transform of Stp(x, t), and 
the total emission probability is 



Ptot(r) = / dk 5<4>{k, t, t -t oo) 



(7) 



We model the target atom based on the soft-core 
Coulomb potential 



U(x) 



1 



\Jx 2 + a 2 



(8) 



and adjust the parameter a = v2 to the ground state 
binding energy eb — —13.6 eV of the hydrogen atom. 
We refer to the exact solution of ©-|g|) as "TDSE re- 
sult" and retrieve the SEA results by ignoring U(x) in ([3]). 
The comparison of TDSE and SEA results is shown in 
Fig.sQJa)-(c). By dropping the laser interaction xE^t) 
in ([2]), we verified numerically that for the given param- 
eters the polarization of this initial state by the ir pulse 
can be neglected [l9[ . The initial state polarization ef- 
fect on the streaked xuv PE spectrum is further discussed 
in Sec. HVl 



III. EIKONAL APPROXIMATION 

In order to trace the influence of the combined action 
of Coulomb potential and ir pulse on the PE, we write 
the PE wave function as 

4> k (x,t) = a k {x,t)e^ k+AL ^ x - ikH l 2+iS ^ (9) 



with a local phase Sk(x,t). The real amplitude a k is 
not important for the present investigation. In SEA, the 
phase Sk(x, t) is given by the Volkov phase 



Sl FA {t) =kj dt'A L (t') 



(10) 



and independent of x. In EA, and without the ir field, the 
phase accumulated by the PE during its propagation in 
U(x) from the location x at time t to the electron detector 
is calculated along the free-electron classical trajectory 
x'(t',t,x) =x + k(t'-t) [H 

/OO i pOO 

dt'U[x'(t , ,t,x)} = - J dx'U(x'). (11) 

In the presence of the ir field, the free-electron classical 
trajectory is modified by a laser-induced drift 

Ax{t,t')=J dt"A L (t") (12) 



to become 



xz,(t , t, x) = x (t , t, x) + Ax(t, t ). 



(13) 



Replacing x' with x^ in (|11|) . we obtain the CL phase [Tl|, 



X,t) 



dt'U[x L (t',t,x)], 



(14) 



and an eikonal approximation to the local phase in ([9]) 

Si A (x,t)=Sl FA (t)+S^ L (x,t). (15) 

In typical streaking experiments and for this study, the 
ir intensity (~ 10 12 W/cm 2 ) is low enough for Ax(t,t') 
being a small deviation from x'(t',t,x). We thus expand 
S£ L about x'(t', t, x) and obtain to first order in Ax(t, t') 

/oo 
dt'F[x'(t',t,x)}Ax{t,t'), (16) 



with the Coulomb force 



F[x'(t',t,x)} 



dU[x'(t',t,x)} 
dx' 



(17) 



The first term in (|16[) is the laser-free eikonal Coulomb 
phase. This term is independent of time as we explicitly 
indicate in (TTT|) . As we should show below this phase 
does not induce any temporal shift in the streaked xuv 
spectrum, but changes the transition probability. The 
second term in (fT6|) includes Coulomb scattering of the 
PE while it absorbs or releases ir photons [l7|- It is pro- 
portional to the ir vector potential and causes a temporal 
shift in the streaked spectrum. 

The numerical results for ([T6")) in Fig. [5] (a) resolve 
the spatial contributions to the CL coupling phase 
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FIG. 2: (Color online) (a) Time evolution of the CL coupling 
phase (see text). Contribution to the PE streaking 5Eqq E (t) 
at x — due to final state CL coupling at (b) hwx = 90 and 
(c) 25 eV. 

S% L (x,t) - S%(x). Keeping in mind that the EA is de- 
signed for short PE de-Broglie wavelengths [l6[ , we con- 
firmed by comparison with full TDSE results (not shown) 
that S^ L (x,t) remains appropriate down to k = 1, i.e., 
a PE energy of « 14 eV This supports the validity 
of the EA for the range of PE kinetic energies in Fig.sQ] 
and [2] 

The transition amplitude for xuv photoemission from 
the initial state ipi to the final state ip k , 



T k {r) = -i dt (r k (t)\xE x (t + r)\Mt)}, 



(18) 



provides the PE probability P(E = k 2 /2,r) = \T k {r)\ 2 
as an alternative to ((5|). Neglecting the laser distortion 
of the initial state (using ipi(x,t) sa ip i (x)e~ l£Bt ) and 
employing the EA-approximated PE wave function for 
4>k{x,t), we obtain 

Tk A (r) =-i J dtdx a(x,t)e-^ k+ALW]x x?P(x) 

xE x (t + T ) e -^ E k A t*t) R ~^ B -^l2)t^ (ig) 

The COE of the spectrum for a free PE would be Eco e = 
fc 2 /2 = ujx — |es|- The (local) energy shift caused by the 
ir field and CL coupling in EA is given by 



J COE 



(x,t) = dS£ A (x,t)/dt, 



(20) 



and does not depend on the time- independent laser-free 
cikonal phase S k (x) in (fT6|) . For sub-fs xuv pulses, con- 
tributions to the time integral (j!9p mainly arise near 
the center of the xuv pulse at t = — r. Approximating 
8Eqq E {x,t) w dS(x,t = —r)/dt, we obtain 



SEE a e {x,t) = -kA L {T)+E°%*(x,T)+E%2*{x,T), 

(21) 



CL,2, 



where 

E^(x,t) = ^A l (t), 



E, 



COE 

CL.2 
COE 



(x,t) 



dx'F{x')A L { ^-r^- - 



(22) 
(23) 



arc the two contributions to the CL shift, E^q E = 



r CL,l „CL,1 
J COE "+" ^COEi 

for hw x = 90 and 25 cV. 



Ecoe + Eco'fi, shown in Fig. s El (b) and (c) at 
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FIG. 3: (Color online) (a) Local temporal shift 5t(x), and 
(b) local oscillation amplitude K(x) induced by the CL inter- 
action in EA. Comparison of the streaked COEs in SFA for 
fajjx = 25 eV with streaking energies (c) in EA at x = and 
(d) from full TDSE calculations. The PE is assumed to move 
to the right (k > 0). 

According to (|22|) and ([23)) . Eqq E (x, t) is proportional 
to 1/k, while SE^ a (t) is proportional to k. Therefore, 
the CL coupling effect decreases for increasing PE kinitic 
energies. As shown in Fig. s El (b) and (c), the cance- 

CL 1 CL 2 

lation between E co ' E and E COE becomes stronger and 
further reduces the CL coupling with increasing k. Note 
that Eqq E (x,t) mainly increases the oscillation ampli- 
tude of the COE in SFA, while E^q' e (x,t) changes the 
oscillation amplitude and induces a phase shift. We can 
thus introduce a local temporal shift St(x) (relative to 
the SFA phase) and a local oscillation amplitude K(x) 
by rewriting (|2Tj) as 



SE, 



EA 
COE 



(x, 



K{x)Ai 



8t(x)] 



(24) 



Fig.s[3](a) and (b) show St(x) and K{x) as functions of x 
for different ujx- St(x) can be positive or negative. The 
actual shift St in the streaking spectrum is obtained by 
spatial integration according to (|19|) . However, as shown 



in Fig. [3] (c) and (d), SE^q E (x = 0,r) agrees well with 
the TDSE result, since the initial wave function ip(x) 
is localized at x = 0. Similarly, we find that the full 
TDSE results for St (solid line in Fig. Q] (d)) and K/k 
(solid line in Fig. Q] (e)) agree well with the EA results 
St{x) (dashed line in Fig. Q] (d)) and K(x)/k (dashed 
line in Fig. Q] (e)) evaluated at x = 0. This justifies 
approximating ASj? A {x,t) = Sj? A (x,t) - S%(x) in (fT5"f 



ASr(x,t)*S^ A (t) 



dt'F[x'(t',t,x = 0)}Ax(t,t') 



KE L {t + 5 T ) 



(25) 



where, in the second line, the slow varying envelop ap- 
proximation is used. 

We notice from (J2TJ) - (J23J) that the three contributions 
to SEgQ E (x, t) are equally proportional to the ir electric 
field amplitude. Therefore, St and K/k do not depend on 
the intensity of the ir pulse. However, reducing ojl, \St\ 
increases and K/k decreases (Fig. [J). This is consistent 
with (|2U|) : at smaller uj^, A^ oscillates slower, leading 
to less cancelation in the time integral and thus to larger 
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FIG. 4: (Color online) Comparison of eikonal temporal shift 
St and oscillation amplitude ratio K/k at two ujl- 
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FIG. 5: (Color online) Polarization effect of the initial state 
on ir-streaked photoemission from the ground state of Id 
model hydrogen atoms for fkdx = 25 eV. (a) Spectrogram 
with initial state polarization, (b) Spectrogram without ini- 
tial state polarization, (c) Corresponding centers of energies 
5Ecoe{t). 

C L 1 

St. Simultaneously, stronger cancelation between SE co ' E 
and 5E C Qg, results in smaller K(x). 



IV. POLARIZATION OF THE INITIAL STATE 
BY THE IR PULSE 

The effect of initial state polarization by the ir pulse 
on the streaked xuv photoemission spectrum has been 
addressed previously |19l - [2l| |. In this section, we analyze 
how it affects the temporal shift St and the oscillation 
amplitude ratio K/k. We find that the significance of 
the initial-state polarization depends on whether or not 
the laser unperturbed initial state is energetically isolated 
from other levels. 

We first consider the non-degenerate case. For our 
one-dimensional hydrogen atom, all levels are non- 
degenerate. In our TDSE calculation, the initial-state 
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FIG. 6: (Color online) Same as Fig. [SI but for the first excited 
state of the Id hydrogen atom. 

polarization by the ir pulse can be included (excluded) 
by keeping (dropping) the term xE^t) in @. In Figs. [5] 
and [HI we compare the polarized (a) and unpolarized (b) 
spectrograms and their corresponding centers of energy 
(c) for the ground state level and the first excited state. 
Due to its large separation in energy from all excited 
states, the effect of the polarization of the ground state 
by the laser pulse on the spectrum is small. It slightly in- 
creases the oscillation amplitude but barely changes the 
temporal shift St. In contrast, the first excited state, 
whose binding energy is 6.34 eV, can be easily polarized 
due to its laser-induced coupling to the second excited 
level at 3.64 eV. As can be seen in Fig. [6j the polariza- 
tion distorts the spectrogram for negative delays where 
the ir pulse precedes the xuv pulse. This distortion does 
not uniformly shift the spectrograms. Therefore, St can- 
not be uniquely defined as a delay-independent temporal 
shift. 

For the degenerate case, we consider a space spanned 
by the stationary wavefunctions ^200 (r) and ^210 (r) of 
the real (3-dimensional) hydrogen atom. Under the in- 
fluence of the laser pulse, the wavefunction is 



ip(r,t) = ^ooW^ooM +0210(^210 (r)] 



(26) 



By shifting the energy scale such that the binding ener 
gies of the two degenerate stationary states are £200 = 
£210 = 0, and by substituting ip(r,t) into the TDSE 

. d 



1— V>(r,t) = [H at + zE L {t)}Hv,t), 



(27) 



we obtain the equations of motion for the coefficients 
a 2 oo(t) and a 2 io(i), 



i Jt a 200 (t) 
l Jt a2w[t) 



nE L (t)a 2 w(t), 
fJ.E L (t)a 2 oo(t), 



(28) 
(29) 



where H at is the atomic Hamiltonian and (i = 
( - 02oo( r )| ;2: IV'2io( r )) the dipole-coupling matrix element. 
The above equations can be solved analytically f22j , 

0200 (t) = a°oo cos [fiA L (t)] + ia Q 2W sin \\iA L (*)] , (30) 
a 2 io(t) = a 210 cos \jj,A L (t)\ + ia,2 00 sin [(J.A L (t)) , (31) 
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where CI2Q0 and a^n are the initial amplitudes at to 



-00. 



For example, the initial values a§ 00 = ±l/\/2 and 
a 2io = l/v2 give the wavefunctions 



ip±(T,t) = V>±(r)e 



which evolve from the Stark states 

^±(r) = [V>2io (r) ± -0200 (r)] 



(32) 



(33) 



Similarly, for a^oo = 1 and cQw = 0, we obtain the wave- 
function 

lMr,t) = -L [v>+(r)e^« - V-(r)e-^W] , (34) 

which evolves from an initial 2s state, while a^oo = and 
a 2io = 1 results in a wavefunction that evolves from an 
stationary 2p state, 

^2p(r,t) = 4= U+We^W +V-(r)c _<Mi(t) l ■ (35) 
\/2 L 
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FIG. 7: (Color online) Photoelectron spectrum for an initial 
2s state of hydrogen: (a) Neglecting ir-laser-induced initial- 
state polarization by setting fi = in (|34|) . (b) Including 
laser-induced initial-state polarization, using fi = 3 a.u. in 
(|34[) . (c) Corresponding centers of energies Ecoe(t). No 
polarization-induced temporal shift is observed between the 
results for polarized and unpolarized initial states. Similar 
results (not shown) were obtained for initial 2p states, using 
(|35l). 



Next, we calculate the ir-streaked spectrum using any 
of the wavefunctions ((32|) - ([35|) as the initial state in (fl8|) . 
In order to disentangle temporal shifts induced by i) the 
initial-state polarization (relative to an unpolarized tar- 
get) and ii) the Coulomb potential acting on the final 
PE state (relative to the SFA, see section III), we ne- 
glect the final-state distortion by the Coulomb poten- 
tial and study initial-state polarization effects within the 
SFA. In this polarization-effect study, we hence use the 
Volkov wavefunction ifa (r, t) as an approximation to the 
final state. If the initial state evolves from a station- 
ary 2s or 2p state according to (|3~4"1) or ([331) . the po- 
larization causes delay- dependent interferences between 
the two Stark states ip±. This interference significantly 
changes the energy-differential PE yield in the streaking 
trace in Fig. [7] b relative to the trace for an unpolar- 
ized initial state in Fig. [7] a. However, the interference 
does not induce a relative temporal shift of the polar- 
ized relative to the unpolarized spectrum, which is best 
seen in the centers of energy of the two spectra in Fig. [7] 
c. This lack of an interference-induced temporal shift is 
explained by the fact that the dipole expectation values 
(4>2s{r,t)\z\ip2s{r,t)) and (tp 2 p(r,t)\z\ip 2p (r,t)) are zero 
at all times, even though the ir-laser pulse mixes the sta- 
tionary 2s and 2p states. 

The situation is different for the states (|32|) that evolve 
out of initial Stark states ip±. The comparison of the 
energy-differential PE yields in the streaking trace for 
initial states (|32|) with and without including ir-laser- 
induced initial-state polarization shows only very small, 
hardly noticeable, differences (Fig.s [5] (a-c)). However, 
temporal shifts [2l| become noticeable in the correspond- 
ing centers of energies. At a PE energy of 60 eV, they 
amount to 41 as between an ir-laser polarized and unpo- 
larized initial ip+ state (Fig.[5](d)) and to 82 as between 



polarized initial %p + and -0_ states (Fig. [5] (e)). These 
shifts originate in the permanent dipole moments of the 
Stark states whose interaction with the ir-laser electric 
field shifts the streaked spectra. 
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FIG. 8: (Color online) Photoelectron spectra for the initial 
n=2 Stark states (|32[) of hydrogen: (a) For 0+(r, t), neglect- 
ing ir-laser-induced initial-state polarization by setting /i = 
in (I32I) . (b) For ip+(r, t), including laser-induced initial-state 
polarization, using /i = 3 a.u. in (|32[) . (c) For ip- (r, t), includ- 
ing laser-induced initial-state polarization, using /j, — 3 a.u. in 
(|32[) . (d,e) Corresponding centers of energies Ecoe(t), show- 
ing a relative temporal shift between the streaking traces of 
(d) polarized and unpolarized initial %))+ (r, t) states and (e) 
polarized initial ^+(r,t) and t[>-(r,t) states. 
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V. CONCLUSIONS 

We have shown how the simultaneous interaction of 
an xuv PE with the electric field of a streaking ir laser 
pulse and the Coulomb potential of the residual ion in- 
duces a specific Coulomb-Laser-coupling phase and leads 
to an attosccond temporal shift and amplitude enhance- 
ment in the oscillation of the streaked PE spectrum. This 
shift and amplitude enhancement become significant and 
observable as the xuv photon energy approaches the ion- 
ization threshold. It can be explained semiclassically in 
terms of an added Coulomb-phase factor in the PE wave 
function. This factor reveals the origin of the observable 
temporal shift as a Coulomb-laser coupling effect in the 
PE dynamics: the PE absorbs and releases ir photons 
while moving subject to the ionic Coulomb force. The 
analytical results obtained in EA show that the CL cou- 
pling induces a temporal shift relative to Al, thus rela- 
tive to the SFA result. For the experimental observation 
of St and K/k as a function of the PE kinetic energy, 
we suggest using xuv pulses with tunable xuv photon en- 
ergy to photoemit electrons from two levels with a large 
energy separation [23j . 

I 



We have also examined the effect of ir-laser-induced 
polarization of the initial state on the ir-streaked xuv PE 
spectrum. If the initial state is not degenerate and has a 
large energetic separation from all other states, its very 
small polarization does not noticeably affect the PE spec- 
trum. On the other hand, if the initial state can easily be 
coupled to other states by the ir-lascr pulse, its polariza- 
tion is important and, interestingly, does not uniformly 
shift the spectrum. If the initial state has a permanent 
dipole moment, such as the n=2 Stark states of hydro- 
gen, there is a relative temporal shift in the streaking 
traces i) for different initial Stark states and ii) with and 
without inclusion of the initial-state polarization. 
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Appendix A: Using Eikonal Wavefunction for RABITT 

In this appendix, we show that atomic phase in RABITT and the relative temporal shift induced by the Coulomb in- 
teraction are identical within an eikonal approximation. We start from the ir assisted single xuv-photon photoemission 
amplitude 

/+oo 
dt{tt>f(t)\T-E x (t + T)\il>i(t)) (Al) 
-oo 

where tpf{t) and tpi(t) are ir-dressed final and initial states, respectively Different from streaking, in RABITT 
attosccond pulse trains (APT) synthesized from a number of odd harmonics of the ir field 

E x (t)=E Xt0 J2e- iu » t+i ^ (A2) 

n 

are used. Here u>„ = (2n + 1)lol and tp n are the frequency and phase of the (2n + l)-th harmonic, respectively, and ojl 
is the fundamental frequency of the ir field. For simplicity, all harmonics are assumed to have the same strength Ex,o- 
Using these APT, we obtain a series of peaks in the PE spectrum that are separated by twice the ir photon energy. 
In the presence of a weak ir field, sidebands will form between the main peaks due to the emission or absorption 
of ir photons. In RABITT periodic intensity variations are observed in the first sideband due to the inteference of 
two distinct two-photon transition routes: 1) absorption of one harmonic photon with frequency u n = (2n + 1)ujl 
and emission of an ir photon with frequency ojl, 2) absorption of an adjacent lower harmonic photon with frequency 
u> n —i = (2n — 1)oJl and an ir photon with frequency u>l. The intensity of the sideband is controlled by the delay r 
between the APT and the ir field, 

P sb (r) ~ [1 - cos (2lu l t + Aif n - A0 Qt )] . (A3) 

It is shifted by the harmonic phases Aip n = ip n -\ — (p„ and the atomic phase Acf> at [T3, EH, [Hi ■ A<f> at is a function of 
the PE energy. 

Using the eikonal approximated wavefunction, we can now show that the atomic phase A(j) at /2ujl is equal to the 
temporal shift St in streaked spectra. Using Eq. (|2"5|) in the main text for the EA phase S^ A , expanding ^ A (x,t) 



up to first order in El, 



(2tt) 1 /2 



e i(k+A L (t))x e -iE k t+iS c (x) 



iKE L (t + Sr) 



(A4) 



substituting it into (|Aip for ipf, and carrying out the time integration, we obtain the transition amplitude up to 
two-photon process 

T(r) - - id k E Xfi ^ S [E + I p - (2n + l)w L ] 

- Z ^^M V {s(E k +I P - 2nu L )e-***-V -6[E k + I p - (2n + 2) Wi ] e^"^} , (A5) 

Wr Z ^ J 



where Ek = fc 2 /2 is the PE kinetic energy, e *(- E fc+-*p) T has been dropped, 



and 



e l{u) = [S(u + u> L ) - S(u> - U> L )] , 
2,1 

E x (uj) = J2 E x,oe +i<Pn 5[u> - (2n+l)w L ] 



(A6) 

(A7) 
(A8) 



are the Fourier transformations of El{£) and Exit)- The transition amplitude for sidebands at PE energies Ek 
2nuiL — Ip follows as 



T 3b (E f ,r) = - i Kd _ kE n L -° Ex.oe-^-^- 6 ^ 



2u>* 



1 - e 



iA(p n -\-2i(jj£ l (t — 6t) 



Accordingly, the sideband intensity as a function of delay r is 



P sb (E k ,r) = \T sb (r)\ 2 = - dk ^ L f Ex >° {i _ C os [2w L (r - St) + Atp n ]} 

ZuJ L 



(A9) 



(A10) 



from which the atomic phase A<j) at /2uil = St is identified. 
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